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Abstract 

It has been shown that finding generic Bell states diagonal entanglement witnesses 
(BDEW) for di ® (^2 ® •••• ® systems exactly reduces to a linear programming if the 
feasible region be a polygon by itself and approximately obtains via linear programming 
if the feasible region is not a polygon. Since solving linear programming for generic case 
is difficult, the multi-qubits, 2 N and 3 ® 3 systems for the special case of generic 
BDEW for some particular choice of their parameters have been considered. In the rest 
of this paper we obtain the optimal non decomposable entanglement witness for 3 (8) 3 
system for some particular choice of its parameters. By proving the optimality of the 
well known reduction map and combining it with the optimal and non-decomposable 
3 (8) 3 BDEW (named critical entanglement witnesses) the family of optimal and non- 
decomposable 3 (g) 3 BDEW have also been obtained. Using the approximately critical 
entanglement witnesses, some 3 3 bound entangled states are so detected. So the well 
known Choi map as a particular case of the positive map in connection with this witness 
via Jamiolkowski isomorphism has been considered which approximately is obtained via 
linear programming. 

Keyw^ords: Entanglement witness, Bell decomposable state, non decom- 
posable entanglement witness, Optimal entanglement witness, Choi map. 
PACs Index: 03.65.Ud 
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Entanglement is one of the most fascinating features of quantum mechanics. As Einstein, 
Podolsky and Rosen [1] pointed out, the quantum states of two physically separated systems 
that interacted in the past can defy our intuitions about the outcome of local measurements. 
Moreover, it has recently been recognized that entanglement is a very important resource in 
quantum information processing [2]. A bipartite mixed state is said to be separable [3] (not 
entangled) if it can be written as a convex combination of pure product states. 

A separability criterion is based on a simple property that can be shown to hold for every 
separable state. If some state does not satisfy this property, then it must be entangled. But 
the converse does not necessarily imply the state to be separable. One of the first and most 
widely used related criterion is the Positive Partial Transpose (PPT) criterion, introduced by 
Peres [4]. Furthermore, the necessary and sufficient condition for separability in H2 <8) H2 and 
H2 (8) was shown by Horodeckis [5] , which was based on a previous work by Woronowicz [6] . 
However, in higher dimensions, there are PPT states that are nonetheless entangled, as was 
first shown in [7], based on [6]. These states are called bound entangled states because they 
have the peculiar property that no entanglement can be distilled from them by local operations 
[8]. 

Another approach to distinguish separable states from entangled states involves the so 
called entanglement witness (EW) [9] . An EW for a given entangled state p is an observable 
W whose expectation value is nonnegative on any separable state, but strictly negative on an 
entangled state p. 

There is a correspondence relating entanglement witnesses to linear positive (but not com- 
pletely positive) maps from the operators on Ha to the operators on Hb via Jamiolkowski 
isomorphism, or vice versa[10]. 

There has been much work on the separability problem, particularly from the Innsbruck- 
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Hannover group, as reviewed in [11, 12], that emphasizes convexity and proceeds by character- 
izing entanglement witnesses in terms of their extreme points, the so-called optimal entangle- 
ment witnesses [13], and PPT entangled states in terms of their extreme points, the edge PPT 
entangled states [14, 15]. 

Having constructed the EW , one can decompose it into a sum of local measurements, 
then the expectation value can be measured with simple method. This decomposition has 
to be optimized in a certain way since we want to use the smallest number of measurements 
possible[18, 19, 20, 21]. 

In this paper, we show that finding generic Bell states diagonal entanglement witnesses 
(BDEW) for di (g) (^2 ® ■■■■ <8) c^n systems reduces to a convex optimization problem. If the 
feasible region for this optimization problem constructs a polygon by itself, the corresponding 
boundary points of the convex hull will minimize exactly optimization problem. This problem 
is called linear programming , and the simplex method is the easiest way of solving it. If the 
feasible region is not a polygon, with the help of tangent planes in this region at points which 
are determined either analytically or numerically one can define new convex hull which is a 
polygon and has encircled the feasible region. The points on the boundary of the polygon can 
approximately determine the minimum value of optimization problem. Thus approximated 
value is obtained via linear programming. In general it is difficult to find this region and solve 
the optimization problem, thus it is difficult to find any generic multipartite EW. In the fol- 
lowing sections we consider some simple but important examples which are solved exactly or 
approximately by linear programming method. Then we consider the multi-qubits and 2<S) N 
with exactly minimum value by linear programming and 3(8)3 systems with approximately min- 
imum value by linear programming and then establish 3 ® 3 optimality condition together with 
non-decomposability properties for some particular choice of its parameters. Then we combine 
the optimal well known reduction map, and the optimal as well as the non-decomposable 3 
(g) 3 BDEW (i.e., the critical entanglement witnesses) to obtain further family of optimal and 
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non-decomposable 3 ® 3 BDEW. Finally, using the critical entanglement witnesses some 3 (8) 3 
bound entangled states are detected and we consider the well known Choi map as a particular 
case of the positive map in connection with this witness via Jamiolkowski isomorphism which 
approximately is obtained via linear programming. 

The paper is organized as follows: 
In section 2 we give a brief review of entanglement witness. In section 3 we show that finding 
generic Bell states diagonal entanglement witnesses for di <^ d2 <^ ■■■■ <S> dn systems reduces to 
a linear programming problem. In section 4, we consider BDEW for multi-qubit system. In 
section 5, we provide BDEW for 2 ^ N. In section 6, we provide BDEW for 3 3 systems. 
Section 7 is devoted to prove the n-d of critical EW and introduce a new family of optimal 
nd-EW via combining critical EW with the well known reduction maps. In section 8, using the 
critical EW, we will be able to detect a bound BD entangled state. In section 9, we consider 
the well known Choi map as a particular case of the positive map connect with this witness via 
Jamiolkowski isomorphism. Finally in section 10 using the optimal EW, we show that some 
separable Bell states diagonal lies at the boundary of separable region. The paper is ended with 
a brief conclusion together with three appendices devoted to the proof of A) the optimization 
of product distributions B)optimality of critical, reduction map C)simplex method for solving 
linear programming problem. 



2 Entanglement witness 

Here we mention briefly those concepts and deflnitions of EW that will be needed in the sequel, 
a more detailed treatment may be found for example in [6, 10, 17]. 

Let S be a convex compact set in a finite dimensional Banach space. Let p be a point in 
the space with p which is not in S. Then there exists a hyperplane[17] that separates p from 
S. 
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A hermitian operator (an observable) W is called an entanglement witness (EW) iff 



3psuch that Tr{pW) < 



(2-1) 



Vp' e S Tr{p'W) > 0. 



(2-2) 



Definition 1: An EW is decomposable iff there exists operators P, Q such that 



W = P + Q^^ P,Q>0. 



(2-3) 



Decomposable EW can not detect PPT entangled states[6]. 

Definition 2: An EW is called non-decomposable entanglement witness (nd-EW) iff there 
exists at least one PPT entangled state which the witness detects [6]. 

Definition 3: The (decomposable) entanglement witness is tangent to S (P) iff there exists 
a (7 e ,5 ( p e P) with Tr(Wa) = (Tr(Wp) = 0). 

Using these definitions we can restate the consequences of the Hahn-Banach theorem [17] 
in several ways: 

Theorem: 

1- p is entangled iff there exists a witness W such that Tr{pW) < 0. 

2- p is a PPT entangled state iff there exists an non-decomposable entanglement witness 
W such that Tr{pW) < 0. 

3- a is separable iff for all EW Tr{Wa) > 0. 

From theoretical point of view this theorem is quite powerful. However, it is not useful to 
construct witnesses that detect a given state p. 

We know that a strong relation was developed between entanglement witnesses and positive 
maps[6, 10]. Notice that an entanglement witness only gives one condition (namely Tr{Wp) < 
0) while for the map (I a (8) 0)^ to be positive definite, there are many conditions that have to 
be satisfied. Thus the map is much stronger, while the witnesses are much weaker in detecting 
entanglement. It is shown that this concept is able to provide a more detailed classification of 
entangled states. 
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As we know, one can expand any trace class observable in the Bell basis as 



W ^ y W. . . lib. . . ){ib. . 



(3-4) 



where iV^i^ij ^ ^ < ^2 < (i2, < i„ < and d\<d2< ... < dn) stands for the 

orthonormal states for a di d2--- <H) (^n Bell state defined as 



where ft and S are phase modules and shift operators for a o?i (8) c?2 <5 



(3-5) 



dn defined as 



n 



1 
a; 








y ... uj"^-^ J 



. S 



1 
1 








y 1 ... y 



(3-6) 



with oj = exp{-j') and 



yd i=o 



(3-7) 



W is a trace one observable i.e., Tr(W) = 1 and we have J] . . . =1. 
Let us split the observable W into its positive and negative spectra as: 

n+ n~ 



(3-8) 



fc=l k=l 

where A^(A^) are the positive (negative) eigenvalues \4>k)i\4>k))j have n'^ + n~ = d^. 

Denoting | |= s > we can write (3-8) as: 



= (l + s)p+-sp-. 



(3-9) 



where p are two normalized positive operators or density matrices defined as 

* fe=l * k=l 



(3-10) 
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Now, using the Lewenstein-Sanpera technique [25, 26, 27] the identity operator ^^^^^ '^^ can 
be written in terms of p~ and the other positive states as 



Idid2...dn 



= Xp- + (1 - X)p' , < A < 1. (3-11) 



did2---dn 

By using the above equation we can replace p" in Eq.(3-9) in terms of the identity operator. 

So, Eq.(3-9) is written as a sum of the identity and positive operators. Thus we have 

W = r^^^ + (1 - r)p, (3-12) 
did2...dn 

where 

[l + s)X , ^ ,1-A 



p+ + s{—^)p-\ (3-13) 



X + s ^s + X 



and r = -J <0. 



In this paper we have considered only trace one observables which are diagonal in the Bell 
states. Hence we restrict ourselves to the Bell states diagonal p defined as 

p= y q . ^ . ){th. . . I , g. , . > and V g, . = 1. (3-14) 

Finally, by substituting (3-14) in (3-12) the trace one Bell states diagonal W observables are 

^^ h,d2...d y q. . .1^,. . )(^. . .1. (3-15) 

^ ^'^ " <li2--in 

The observable given by (3-15) is not a positive operator and can not be an EW provided 
that its expectation value on any pure product state is positive. For a given product state 
It) = |q;)i|q;)2...|q;)„ the non negativity of 

Tr{W\^){^\) > (3-16) 

implies that 

- <r < 0, (3-17) 



1 - did2...d„T,. . . q P 
where P . =|< 7 1'?/' . >P. 
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Denoting the summation in the numerator and the dominator in (3-17) by ^"(7) — 

dido-.-dnY,. . ■ Q P we see that the least possible rn = — t^t^tt is the decreasing 
function of C(7) for C(7) < 1 (obviously for C(7) > 1 all r while being positive provide positive 
expectation value). Therefore, for given parameters q. . . > ,with y;. . a = 
1 , the least allowed value of the parameter r, called the critical parameter (denoted by 
Yc ) is obtained from the product state 7 which minimizes C-v = 1]. . - Q P . , with 
< P. . . < 1 and the constraint y^. . . P. . . = 1. As for the completeness of the 
Bell state yi. . IV" 1 = 1, the determination of Tc reduces to the following 

optimization problem[24] 

minimize C-v = y (? P (7) 

< P. . . (7) < i- (3-18) 

^. . . P. . . (7) = 1. 

Always the distribution P.^.^ satisfies < P.^ (7) < ^ for all pure product states (the 
proof is given in the Appendix A). One can calculate the distributions P.^.^ i„(7)) consistent 
with the aforementioned optimization problem, from the information about the boundary of 
feasible region. To achieve the feasible region we obtain the extreme points corresponding to 
the product distributions Pji2...i„(7) every given product states by applying the special 
conditions on q.^-^ -Js parameters. themselves are functions of the product distributions, 
and they are in turn are functions of 7. They are not real variables of 7 but the product states 
will be multiplicative. If this feasible region constructs a polygon by itself, the corresponding 
boundary points of the convex hull will minimize exactly in Eq. (3-18). This problem is 
called linear programming , and the simplex method is the easiest way of solving it. If the 
feasible region is not a polygon, with the help of tangent planes in this region at points which 
are determined either analytically or numerically one can define new convex hull which is a 
polygon and has encircled the feasible region. The points on the boundary of the polygon can 
approximately determine the minimum value C-y from Eq.(3-18), thus the problem is that of a 
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linear programming again. In general it is difficult to find this region and solve the optimization 
problem, thus it is difficult to find any generic multipartite EW. In the following sections we 
consider some simple but important examples which are solved as linear programming problem. 



4 Bell states diagonal entanglement witnesses for multi- 
qubit system 

Here we provide a multi-qubit entanglement witness. From the previous section one can show 
that the Bell states diagonal observable W for multi qubit system is defined by 

I n ^ 

W = r^ + (l-r) Yl fc,«2,...,inlV'n,i2,...,in)(V'n,i2,...,i„l, (4-19) 
n,...,in=0 

where |'?/'ii,i2,...,i„) is a Bell state: 

= ® ® - ® M''\^o,o,...,o): (4-20) 

with 

Ko,...,o) = ^EK)lK)2-K)n, (4-21) 

and and are the Pauli operators. This observable is not a positive operator and can not 
be an EW provided that its expectation value on any product state I7) = |q;)i|q;)2...|q;)„ is 
positive. 

We consider an easy case q^^ — , q^^ — x with all the other q^s being equal, i.e., 

^ii,»2--.«n = 2(2"- ^"-1) except for = = ... = in = and ^2 = is = ■■■ = in = 0, ii = 1. Then 
the observable W reduces to the following form 

^ = r ^ + ^(^^r^( (1 - ^)^2" - (1 - ^) l^o,o,...,o) (^^^^^^^^ 

(4-22) 

We can calculate from the non negativity of T'r(l^|7)(7|) for a given product state I7) 

= 2(2n-^i_l) ((^ - ^) - (1 - ^)^oo...oo + ((2" - 1)^ - mo...oo)- (4-23) 



Bell states diagonal entanglement witnesses 



11 



According to the definition of product distributions, we liave 



(4-24) 



where 



\a) 



Vft; 



, i = l,2,...,n. 



(4-25) 



For the most general given product states 7, we determine the extreme allowed values of these 
product distributions. Let 



\a) 



then the product distributions are 



P. 



_ 1 
00. ..00 2 



P = i 
^ 10. ..00 2 



Another choice will make P^g ^ minimal 



and the product distributions will become 



1 

V2 



(4-26) 



(4-27) 



(4-28) 



10. ..0 2" 



P =0 

10... 00 



Similarly we can also make P^q q minimal 



\0i)2 = = ... = \a)n = ^ 



^ 1 ^ 



p =0 

00. ..00 



1 



10. ..00 2"-l 



(4-29) 



(4-30) 



Prom the definition of the convex function [24] we can show that the convex combination of 
these distributions provide a convex region called the feasible region, where all points in the 
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interior of this region satisfy the positivity constraint of Tr( iy|7)(7|). Then we have 

maximize -C^ = 2(2^-1,1) (-(1 - x) + {1 - x)P^^^^^^^ - ((2" - l)x - l)P,o,„J 
subject to 2Poo 00 - 2Pio 00(1 - 2^^) < 

2P — 2P (^ ^^'i < ^ 

•^-^ 10...OO OC.OOV-^ 2"-!/ — 2"-l 

P > P > 

00. ..00 — '^1 ^ 10. ..00 — 

Now we must prove that the feasible region, constructed from the convex of these points, 
is a polygon. Let P+ = -P00...0 P- — Pio...o- The equation of the hne passing through 

(P+ = ^,P_ = 0) and (P+ = iP- = i) is 

on— 1 1 

P_ = (—- )P+ ^ . (4-32) 

Let us further assume P_ = AP+. By intersecting this equation with the one above we get 

_ 1 
^+ " 2n-i_A(2-i-2)- 

Now if we assume A = we arrive at the point (P+ = P_ = 0), for A = 1 we conclude 

(P+ = |,P_ = I). One can rewrite Eq.(4-24) as 
1 



P± = T^il "1 I I "2 I ••• I a-n I +(1- I Oil I )(1- I Oi2 I )•••(!- I a 



±2 I cci I (1- I «! 1^) I 0L2 I (1- I 0L2 P)... I q;„ I (1- I OLn f) cos(0)). (4-34) 
Thus we write the Lagrangian as 

£ = P+ + /x(P_-AP+), (4-35) 
where \i is the Lagrange multiplier. With | |= cosQi we maximize L with respect to 6'j's and 

01 — 92 — ■■■ — On =^| CKi 1 = 1 a2 \— ... —\ OLn I 



(4-36) 







such that 



tan"^i = ^ (4-37) 
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so that 

P+ = ( — ^) ^ ^ 2 ■ (4-38) 

As we see the equation for P_,. is less than the one in (4-33), moreover, this relation indicates the 
correctness of the result (4-26)- (4-30). Thus the convex hull is a polygon and the optimization 
problem will be converted into the linear programming one. 

There is no simple analytical formula for solving a linear programming, however there are a 
variety of very effective methods, including the simplex method to solve them. So, minimization 
solutions of C-y is obtained by the simplex method[24] and we have (see Appendix C): 

I) For < X < 2n-i+i the extreme points of the feasible region are i-*oo...oo — -P10...00 — i 
the minimum value of C-y is defined by {C^)min = f • By substituting these values in (3-17) we 
have 

— < r < =^ re = 4-39 

where Vc is called the critical r. By substituting Tc in (4-19) this observable has positive 
expectation value under any product state, thus it will be an EW called critical EW equal to 

1 1 — a: (2" — l)x — 1 

Wc{x) = 2(2n-l _ 1) (^2" ~ ^_2n-l^ l^o,o,...,o)(V'o,o,-,ol + ^ _ ^n-l^ IV'i.o.-.o) (V'l.c-.ol))^ 

(4-40) 

which in the special case where x — -^^—^ the Wc{x) reduces to 

^red = 2(2"-^l-l) ^^^" " 2|V'o,o,...,o)(V'o,o,...,ol), (4-41) 

which is the well known reduction map. 

II) For 2^-^+1 ^ ^ ^ 1 the extreme points of the feasible region are P^^ = and 
-^10... 00 — respectively. Therefore, from the simplex method we get {C^)min — hence 
Tc = — i=£ and the critical EW is calculated to be 

/ N 1 ,1 - ,T 1 - a;, , . , , , (2"' - l)x - 1 1 , , , , 

^^(^) = (2„-l _ l) ^~^2^ " ^;^IV'o,o,...,o)(V'o,o,...,ol + ^ \A,o,-,o)Wi,o,...,o\))- 

(4-42) 
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Note that this choice of q is not the only way of defining a BDEW for multi-qubit system in the 
one parameter representation. Let us consider the alternative definition for the one parameter 
BDEW by studying the following example. Assume goo . oi — ^ other q^s to be 

equal. Thus we have 

(4-43) 

Similarly we can find the extreme points of -Poo.. .00 -P00...01 ^ 

= |q;)2 = ■■■ = \oi)n 



P = i 
00. ..00 2 



(4-44) 



P =0 
00. ..01 



\a)i = \a)2 = ■■■ = \a)n-i 



(A 



P =0 
00. ..00 



(4-45) 



and la), 

1 \ I _ 

00. ..01 2 

Also we know that the convex combination of -P00...00 -P00...01 provides a convex region 
or a feasible region. Then we have an optimization problem as follows: 

minimize = 2(2n-i_i) ((l - x) - {I - x)P^^^^^^^ + ((2" - l)x - l)Poo...oi) 



subject to i - - > ■ (4-46) 

P P > 
00. ..00 ' 00. ..01 — " 

Here the optimization is converted to the linear programming problem. To prove, we must 

show that the feasible region is a polygon. Let us suppose P+ = -^00... 00 — -^00... 01 with 



P_ = i I Pia2... an + aiP2---Pn 



(4-47) 



where as before we have P_ = AP+ and | CKj |= cos 9i. By maximizing the Lagrangian we get 

6*2 = 6*3 = • 



9n- Thus we have 



P+ + P- 



1 

2' 



:rlCOs2"-2^2+sin2"-2^2) < ^• 



1 

2" 



(4-48) 
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The line passing through the points (P+ — ^,P- = 0) and (P+ — 0,P- = |) is P_ = —P+ + |, 
which is always located above the curve obtained in the Eq.(4-48). Therefore, all the points are 
within the feasible region and this region constructs a polygon. Thus the above optimization 
problem reduces to a linear programming one. This minimization is exactly solved in the same 
way as mentioned above, and the critical EW is obtained as 

-2"(1 - x) 



2(2"-i(a;+ 1) -2) 



(4-49) 



(2"-i(x + l) -2) 



((1 - a;)J2n - 2(1 - a;)|^„,,„)(^„,,J + 2((2" - l)x - l)|^oo..oi)(^oo..oil)- 



5 Bell states diagonal entanglement witnesses for 20 N 
system 

Here, wc will find a 2 iV entanglement witness. From the previous discussions we can define 
the Bell states diagonal observable W as 

N-l 1 

EZ 

where \i;J ^h® (^)^(Q)°|V'oo), with jV'oo) = ^ELoI^)I^) and 



(5-50) 



10 ••• 

0-1 ••• 

1 ••• 

••• 1 



s 



1 

1 
1 





(5-51) 



/ 



y u u u • • • 1 y y 

Similar to multi-qubit let — and q^^^ — x and let all the other q's be equal to 2n-2- 
Then by obtaining the expectation value of W on the product states and finding the product 



Bell states diagonal entanglement witnesses 



16 



distributions we have 



1 

V2 



( 1 ^ 



F/2 = 7^ 



1 



1 
72 



1 





1 





^0 = 



10 2 



(5-52) 



(5-53) 



(5-54) 



The feasible region is a rectangular and the optimization problem reduces to linear program- 
ming. Therefore, by using simplex method for < a; < we find the minimum value of 
Cj — ^ and critical r as Tc = iz^, and the critical EW is defined as 

1 



2N 



^ ^|y^oo)(V'ool+ ^^^. ^Ir V io)(V'iol)- (5-55) 



For the critical r we find Tc = iii the region < x < 1 and the critical EW has the 

following form 

^'^^^^ = 2^^^ " ^lV^oo)(V^ool + ^ ^ \Ao){^J)- (5-56) 



In another one parameter EW example we assume that q^-^ = x and set all the other g's to 
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be equal so that we have 

Similarly we can find the extreme points of P„„ and as 



(5-57) 



\a)i 



Then the critical EW is defined as 



^0^ 



1 



-2Ar(l -x) 



00 2 



= 



1 



(5-58) 



(5-59) 



{N{x + 1) - 2) 



2(A^(x + l)-2) 

((1 - X)I2N - (1 - ^)|V^oo)(V^ool + ((2A^ - 1)^ - l)lV^oi)(V'o 



(5-60) 



6 Bell states diagonal entanglement witnesses for 3 (8) 3 
system 



Here we provide a 3(8)3 Bell diagonal entanglement witness. One can show that the Eq. (3-15) 
for a 3 (8) 3 system reads as 



I ^ 



(6-61) 



11,12=0 
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It is difficult to prove whether or not the EW for a 3 (8) 3 system is optimaL Also it is difficult 
to see for which value of the allowed r, EW are (or are not) decomposable. Therefore, to 
investigate the optimahty and non decomposabihty of these EW we restrict ourselves below 
to some particular choice of Qij: 

Because the distributions < Pij < | and the minimum value of are dependent on the 
coefficients Qij, we consider a special case for the coefficients qij defined by 

1 1 1 , , 

qoi = qQ2 = qu = q22 = qi2 = q2i = - , qio = x and g2o = ^ - a; , < x < -. (6-62) 

By substituting these values in (3-15) we get 

T 7" 1 Q™ 1 

W{x) = r-| + (1 - r)(^ - -|V'oo)(^ool - -^i\Ao){Ao\ - \Ao){Ao\))- (6-63) 



By using (3-16) for non-negativity of the observable W we find the distributions Pij as a 
function of x. The minimum value of is obtained from the boundary of the feasible region, 
i.e., we have 

C7 = ^(1 - ^0 - (8a; - l)(P,o - ^.o))- (6-64) 

For given product states I7) = |q;)i|q;)2 one can obtain the extreme points of the product 
distributions as 



1 
V3 



\a)i = \a)2 



1 



UJ 



[Pqo — \, Pw — 0, P20 — 0) 



(-P00 — 0, Pio — \, P20 — 0) 



{Poo — 0, Pio — 0, P20 — \) 
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1 

V2 



1 

V2 



, |a)2 = 



1 

x/2 



1 

V2 



^ ^ 



' 1 ^ 



1 1 ^ 



-00 



\ » / 



(Poo — 0, Pio — |, -P20 — \) 



(-P00 — -Pio — 0, P20 — \) 



(6-65) 



(-P00 — -Pio — 75 -P20 — 0) 



and 



3 _ip _ip 

00 — 3) -< 10 — 3) -'20 — 3, 



(6-66) 



By convex combination of these points we obtain the possible region. Thus we have an opti- 
mization problem as 



minimize = \(l - P^^ - {Sx - l)(P,o - P^o)) 
subject to 1 - 3P„„ - P,, + P,o > 

1 + P„„-3P,„-P,„>0 (6-67) 
1-P„„ + P,„-3P,„>0 

P P P > 

00 ' 10 ' 20 — 

One can prove analytically that the region can be encircled with a polygon and the optimization 
problem is reduced to a linear programming. To prove, we begin from the definition of the 
product distributions 

-Poo = I I aiA + (^2^2 + Q;3/33 P 
Pio = I I + a2l32UJ + a^l3^uj p (6-68) 
-f2o = I I Oii(5i a2(520j + as/^so; p . 
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Without loss of generality, one can assume that 

I cti(3i \ = xi , I q;2/92 \ = X2 , I asPa \ = X3. (6-69) 
The Schwartz inequality yields 

Xi + X2 + Xs<l. (6-70) 

Since we are looking the extreme points we will choose the maximum value in the inequality 
(6-70), that is xi -|- X2 + = 1. Thus the product distributions are written as 

Poo^ll xi + X2e''t'' + xse''!'^ p 
p^o = i I + X2e''^^uj + xae'^^oo p (6-71) 
P20 = I I xi X2e''^^LU + x^e^'^^uj p . 

Now, supposing that Pqo ^-nd Pio are fixed values, we conclude 

(p2 = h , X2 = x-i. (6-72) 

Thus 

P00 = I I 1 + 2X2(1 -COS 02) P 

Pio = i I 1 + 2x2(1 - COS02 + f ) P (6-73) 

P20 = I I 1 + 2X2(1 - COS 02 - f ) P ■ 

We write down the equation for the planes passing through the obtained extreme points for 
Pij and maximize it with respect to the variables 02 and X2. Where the obtained values for 
02 and X2 are indicative of the violation from the equations of planes, that is there are points 
which are located out of these planes, or in other word the planes have become convex. Thus, 
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the equations of the planes (Fig 1) and their maximum violation(D) are obtained as follows 

1) 3Pio - P20 + ^'OO - 1 = X2 = ^,COS(/.2 = ^ D=^ 

2) 3Pio + P20 - i'oo - 1 = a;2 = ^,cos02 = ^ = ^ 



3)3P20 + PlO - Poo - 1 = X2 = ^, COS02 = ^ 



4)3P20 - Pio + i^oo - 1 = X2 



61 • 

J_ 

61 



, COS 02 



5)3Poo + P20 - Pio - 1 = cos 4>2 



-1 

7 

13 
14 

13 



P) 

P 

D 



2_ 
61 

_2_ 

61 

_2_ 
61 



(6-74) 



P = ^ 

61 ■ 



6)3Poo - P20 + PlO - 1 = X2^ ^, cos 02 = 14 

It is seen that the points thus obtained are located out of the considered plane. Thus, the 
equations of the planes passing through the new extreme points which are parallel to the above 
plane are obtained. For example, the equation of plane parallel to 3poo + -Pio + p20 = 1 is 
3poo + -Pio + p20 = 1 + ^ which under permutation (Pqo, PiO) -P20) will act similarly. Tack 
arbitrary any three of planes passing through the new extreme, (Pqo = |, Pio = |, -P20 = |), 
(Pio = 0, P20 — 0, Pqo = 0) and Pio + P20 + -Poo — 1 = and intersecting with each other. Hence 
new extreme points will be produced. Thus we have encircled a polygon by its feasible region 
and the optimization problem will be reduced to that of a linear programming. The vertices of 
this polygon are the solutions of the problem provided that they obey (Pqo < |, Pio < |, P20 < 
|) and by substituting them into equation Cj one can determine its minimum value. Thus, 
for ^ < X < the extreme point is defined by (|, |, |) and finally {Cj)min — (^)- Having 
found the critical r we substitute it in (3-15) and obtain a family of EW (called critical EW). 
Thus we have 



Pe = -3 , W,{X) 



-(- 

2^3 



^9 - l^oo)(^ool + (82; - l)(|^io)(V'iol - \Ao){Ao\)), 



(6-75) 



where Wc{x) reduces to the following well known reduced EW at x = |: 

Ig - 31-000) (-000 1 



Wred = 



6 



(6-76) 



In the Appendix B it is shown that the above EW is optimal in contrast to the conclusion that 
it is a decomposable EW Ref . [28] . 
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In the Appendix B, we discuss the possible choice of x consistent with Cmn — ^ and the 
optimahty of the corresponding Wc{x). 

Also we prove in the following section that Wc is nd-EW for all values of < a; < 
except foT X — |. Besides taking a convex combination of Wc and Wred 

WA^AWc+{l-A)Wred, (6-77) 

we obtain a new EW which is optimal (see Appendix B) and is also an nd-EW for certain 
value of the parameter A as will be shown in section 7. 

However we can consider other values for q^j in (3-15), e.g., g2o = Q'02 = Qii = 922 = = 
Q21 — ^ , Qio — X and qoi — \ — x , < x < ^ then define the observable W by substituting 
the above condition in (3-15) as follows 

W{x) = r| + (1 - r)(:| - ll^Ji^J - ^^{\Ao){i^J - (6-78) 

By using (3-16) for non-negativity of the observable W we find the distributions Pij as a 
function of x. The minimum value of is obtained from the boundary of the feasible region, 
i.e., we have 

C, = i(l - P,, - {8x - 1)(P,, - PJ). (6-79) 
For given product states I7) = |q;)i|q;)2 one can obtain the extreme points of the product 
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distributions as 



00 



q;>i 



V3 



V3 



/ - \ 

LO 



UJ 



ah = 



1 

V3 



UJ 



v^y 



a; 



v^y 



v^y 



(-Poo — k, Pio — k, Pqi — 0) 



I! -Poo — k, Pio — 0, -Pqi — I, 



[Poo — 0, -Pio — \, -Pqi — \) 



[Poo — \, Pio — 0, Poi — 0) 



(-Poo — 0, Pio — |, -Poi — 0) 



(Poo — 0, Pio — 0, Poi — h) 



(6-80) 



Similar to (Poo, Pio, -P20) we obtain the extreme points. Then by convex combination of these 
points we obtain the feasible region (see Fig 2) . The optimization problem is in the following 
form 

minimize = |(1 - P^„ - (8a; - 1){P,^ - P^J 



subject to 



~ Poo ~ Pio ~ Poi ^ 



(6-81) 



3 - 00 ~ -T 10 ~ -T 01 ^ 

P P P > 

00 ' 10 ' 01 — 

We have been able to fined analytically the extreme points which at the same time don't 
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violate the plane (| — Pqo — -Pio — -Poi = 0). But we have failed to prove in general that no 

point lies out of the plane. Therefore, we have proved numerically that there is no violation 

from the plane. Thus, this feasible region is a convex hull or a polygon itself and reduces the 

optimization problem to a linear programming. So the vertices of the polygon are the solutions 

of the problem which by substituting them into equation C-y one can determine its minimum 

vahic as (C^)mm = We can find the critical r and by substituting the critical r in 

(3-15) we obtain a family of EW (called critical EW) resulting in 

_ -3 + 24x 
~ -1 + 24x 

^'^""^ ^ ^(ZT^^((8^- l))^9-3|^oo)(^ool -3(8a;- 1)(|^,„)(^,J - |^oi)(^oil))- (6-82) 

The obtained EW for two sets of triplet P, namely (Pqo, -P20, -P20) and (Pqo, -P20, -Poi) can 
produce the most general form EW corresponding to combination of another triplet P^j's. Since 
under Fourier transform one can transform all the shifts in to the modulation. Moreover, the 
shift and modulation operators themselves can affect on the corresponding EW too, and so 
produce new combination of triplet. 

7 Non-decomposible 3 (8) 3 Bell states diagonal entan- 
glement witnesses 

By calculating the partial transpose of PVc(^ ^ x < ^) (for {Poo, -PiO) -P20 case) we prove 
that it is an nd-EW. The necessary and sufficient condition for non-decomposibility of Wc 
reduces to the negativity of its partial transpose. Using the following relation 

(I Vj'fe' >< Ipjk 1)^-^ = ^ IZ^""'! i^m+j',l+k' >< 1pm+j,3-{l-k) I, (7-83) 

l,m 

one can show that (Wc)^^ is a block diagonal, i.e., we have 

{Wcf^ = {Oj)kk' I i^j'k' >< ^jk I, 

j,k,k' 
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with the matrices Oj calculated as 



^ ^ 



i c, 



(7-84) 



with 



(7-85) 



Using the fact that | C2 | = | Ci \ = \ Cq \ one can show that the matrices Oj have the same 
eigenvalues 



A = 



A± = Ai = I ± I J4 + 48a;(4a; - 1^ 



(7-86) 



The above equation indicates that A_ is negative except for the particular case in which a; = |, 
i.e., Wred- Then different eigenvectors are so obtained 



A = 



>=l i^jo >, 



A = A± ^1 (Pf >= iPi I > + I >), 



(7-87) 



where 



± 



Xl-X+B' 



(7-88) 

\± — A±J 

So we conclude that Wj^ has three eigenvalues, namely Aq, A±, each with degeneracy 3, and 
the following projection operators 



Q+ = E-=o I 0/ X 



(7-89) 



Q- = E -=0 1 4>J >< <PJ I 
^ Qo = E|=o I 0° X 0° I • 

Here we have 

W^^ = X+Q+- I A_ I g_. (7-90) 

The equation indicates that W^^ is not a positive definite operator except for the particular 
case Wred: heuce it is non-decomposable entanglement witness. 
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We are interested in the n-d of EW given in (3-15) for the allowed values of p. Therefore, 
we write Eq.(3-15) as 

W ^ sIg/9 + {1 - e)Wc, (7-91) 

with 

. = (T-92) 
Now, expanding Ig/9 in terms of the projection operator (7-89) as 

h^Qo^ + Q-^+Ql^: (7-93) 
the EW given by (3-15) can be written as 

W = e/9Q^' + (| + (1 - ^)^+)^+" + (| - (1 - ^) I ^- 1)^-"- (7-94) 

The above form of EW indicates that its partial transpose W'^^ is positive, i.e., it is decom- 
posable EW if we have 

for ~^^!^"^"| ^ < r < —3. It is not easy to tell where the EW is or is not decomposable. In the 
next section using some bound entangled state we will investigate their non-decomposability. 

Now, in the remaining part of this section we try to obtain some nd-EW by taking the 
convex combination Wc{x) for all ^ < x < and Wred (6-76) as 

Wa{x) = AW,{x) + (1 - A)Wred , A e [0, 1]. (7-96) 

In order to test the positivity of Wj'* {x) we must first expand Wc and Wred in terms of the 
positive diagonal operators. Thus at first we write the projection operators defined in (7-89) 
in the following form 

2 2 

Q± = T.\Xk){Xk\ , Qo^T.\^ko){^ko\ (7-97) 

k=0 k=0 
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with 

1 Xk >= (I ^ki > ±c^~' I ^k2 >). (7-98) 
Now writing Ig/9 in terms of the projection operator (7-97) and using the fact that 

2 2 2 

3 



(I tPoO >< i^OO \ = ^(E I ^kO >< ^fcO I + E I Xfe >< Xfc I - E I Xfc X Xk I) 



fe=0 fe=0 fe=0 

and 

Wj-ix) = A+ E I Xt X Xt I - I A_ I E I Xk X Xk I (7-99) 

k=0 k=0 

we get for the partial transpose WA{x)m Eq.(7-96) 

Wa^-{x) = A(A+) E I xt X xt I +(-A I A- I +^^^) I Xk X Xk I ■ (7-100) 

k=o 

This expression imphes that Wa^^{x) is positive, since 

A < ^7^. (7-101) 

~ 1 + 3 I A_ I ^ ^ 

Again, for | < A < 1, it is not easy to talk about decomposable or non-decomposable 

W\{x), and one needs to find some bound entangled states to show their non-decomposability, 

this will be done in the following section. 



8 Detection of bound entangled state with Bell states 
diagonal entanglement witnesses 

Now if we succeed to find any bound entangled state [6, 5] so that BDEW is able to detect this 
bound state corresponding to BDEW, from definition 2 in section 1 EW will be an nd-EW. 
Let a bound entangled Bell decomposable state be written as 

P = pQo^ + r]Ql^ + CQ-^ , p'^^ > ^ T), C) > 0. (8-102) 

Optimal BDEW must detect this bound state, i.e., 

Tr[W^p\ < ^ 7;A+ < CI A_ I . (8-103) 
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On the other hand this bound state must be positive. For simphcity we use the operator Wc 
and the identity operator /g in the bound state definition 

_ I A- I (79 - Q^^) -W, + XM-Q^-) 

^+ " I A_ I +A+ ' - |;^_ I ' 

so that the bound state reduces to the following form 

|A_|+A, |A_|+A, ^'^+^^rFA;^^- ('-'^'^ 

In this case Qo =| V'oo >< V'oo | + | V^io >< V'lo | + | ^^20 >< ^^20 | and by substituting this 
result in the Eq.(8-105) we get 

V — C V — C 

P={l^- 3(|;)^ I+A+) '' ' ' ^^'^ ^ ^"^^^ ~ ^^ 3(1 A I+A+) '' ' ' 

+ (/. - (12a: - 2)—^^^^) | i^^o >< i^20 |) + ^~ '.T!^ T ^^1 ^^^ (1 V'oi >< V'oi | 



3(|A_|+A+)' 

+ I "002 >< V'02 I + I V'll >< V'll I + I V'22 >< V'22 I + I "012 >< 1pl2 \ + \ 1p21 >< "021 |)- 

(8-106) 

The positivity of p requires that all the Bell states diagonal operator coefficients to be positive, 
and that this condition be imposed on the coefficient fj, only. So we get 



(12a;-l)+3(|A_|+A+) 



(8-107) 

1 (2-12x)(|-2^) 
- 8 ^ - (2-12a;)+3(|A_|+A+)' 

which, in this case means Qo is on the boundary. Now by using this bound entangled BD state 
we can find n-d condition for BDEW. We know EW will be an nd-EW if this EW is able to 
detect any bound state. Then by using the equations (7-94) and (8-102) we have 

Tr{Wp) = {f + 3(^ + (1 - e))rjX^ + 3(^ - (1 - s))C \ A_ |< 0. (8-108) 

Now by substituting e from Eq.(7-92) we obtain 

1 + 27(C I A_ I -^A+) ' 

where the calculated r is greater than the represented r for EW in Eq.(7-95). Therefore, we can 
find one of the p's corresponding to EW which is an nd-EW. Non-decomposable generahzed 
EW for a general case is under investigation. 
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9 Choi map 

Choi positive map [16] 4>{a, b, c) : is defined as 

^ apu + bp22 + cp33 

'^a,b,c{P) = ap22 + bp33 + cpn 

ap33 + bpii + CP22 

where p e M^. It was shown that 4>{a, b, c) is positive iff 

a>l , a + b + c>3 , l<a<3. 



V 



p, (9-110) 



(9-111) 



Using Jamiolkowski [10] isomorphism between the positive map and the operators we obtain 
the following 3 (8) 3 EW corresponding to Choi map 



Wchoi = 



1 



3(a + 6 + c- 1) 



l^fco)('0fco| +bJ2\'^k2){ll^k2\ +cJ2 \i'kl){'^kl\ - 31-000) (-0001 



k=0 



k=0 



k=0 



(9-112) 

Similar to BDEW we expand |'0oo)('0oo| using the identity operator and the other Bell diagonal 
states: 



|^00)(V'00| = ^9 - J2 

Then we reduce EW to the following form 

1 2 

Wchoi = ^, Tt(-(3 - a)h + 3 X! IV'fco)(V'feo| 

6[a-\-o-\-c i.) ^^-^ 

2 2 

+ {b + 3-a)J2 \M{^k2\ + {c + 3-a)J2 \i^ki){Ai\)- 

k=0 k=0 



(9-113) 



(9-114) 



Comparing with BDEW (3-15) we have 



3(3 - a) 



and the EW operator is defined as 



Wchoi^rIg/9 + {l-r){ 



{8-2a + b + c) 



{a + b + c-iy 



l^feo)('0fco| + 



(9-115) 



(6 + 3 -a) 



3(8 - 2a + 6 + c) 



\'^k2){A2\ 
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(c + 3 — a) 



3(8 - 2a + 6 + c] 



J2 IMiM)- 



(9-116) 



fe=0 



By comparing (9-114) with (3-15) we obtain the coefficients qij 

1 b+3-a 



Qio — Q20 



{8-2a + b + c) ' 

qoi = qu = q2i 



qo2 — qi2 — q22 



c + 3 — a 



3(8 - 2a + 6 + c) ' 



(9-117) 



3(8-2a + 6 + c)' 

Note that if r is negative, as introduced in EW above, this operator will be positive, but not a 
completely positive map. For r < we have 1 < a < 3. By assuming a >b > c, the minimum 
negative eigenvalue of choi EW (9-116) is given by 



- + (1- r) + 

9 ^ ^^3(8-2a + 6 + c) 



<0, 



(9-118) 



where upon substituting r from Eq.(9-117) we get 1 < a < 2. This is equal to the introduced 
positivity condition of Choi map in [16]. 

By using (3-16) for non-negativity of the observable Wchoi we find the distributions Pij as a 
function of qij. The minimum value of C-y is obtained from the boundary of the feasible region, 
i.e., we have 

1 _ (6 + 3 -a) _ (c + 3 -a) 



(C7 



_ (6 + 3 - a) _ 



(9-119) 



(8-2a + 6 + c) 3(8-2a + 6 + c) 3(8-2a + 6 + c)' 
where Vi = Z)fc=i -P^o ' ^2 = J2k=oPk2 ^^'^ ^3 — J2k=o^ki- We can find the extreme value of 
{Vi,V2,V3) which is obtained under the product states I7) = |q;)i|q;)2 as 



^ " Oil II A I + I Qi2 II /?2 I e^^2+ I II /?3 I e^^^ |2 

=1 ai n /3i P + I as n /3i p + I as \'\ (3^ \' 
=|ai \'\/32\' + \a2 \'\f3s\' + \as \'\ f3i \' 



3 I 
V2 

Vs 



(9-120) 



where \a)i — 



( \ 



a2 
yas J 



and \a)i — 



One can obtain the extreme points of the 
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(^1,^2,^3) as 



\a)i = \a)2 



la)] 



ah 



a>i 



^0^ 



^ ^ 



(Pl=0,P2 = l,p3 = 0) 



(9-121) 



(P, =0,P2 = 0,-P3 = 1) 



The convex combination of all extreme points provide a convex or a feasible region (Fig-3), 
then we have the following optimization problem 

minimize (C^) = ( (8_2a+fe+e) ^l + 3{£STfc)^2 + 3(iSg^P3 

subject to 1 - fPi - - < 

1 - §7^1 - IP2 - 7^3 < (9-122) 

1 - Pi - P2 - 7^3 > 
Vl,V2,V3 > 0. 

Whether analytically we have been able to show that we will have violation only from the two 

planes 

2 - 3Pi - 2V2 - P3 = 

2 - 3Pi - 2V3 - P2 = 0. 

Now let us assume that the maximum value of the violation from the planes is A < 1. Thus, 
the equation of the plane passing through the new extreme points, parallel to the above plane. 
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is obtained. Next we derive the intersection of the following adjacent planes 

1) 3Vi + P2 + 2P3 - (2 + A) = 

2) 3:Pi + 2P2 + P3-(2 + A) = 
7^1 + 7^2 + 7^3 - 1 = 

Pi =0 
7^2 = 
7^3 = 



3) 
4) 
5) 
6) 
7) 
8) 
9) 



(9-123) 



7^1 = - 

'^i 3 

7^2 = 1 

7^3 = 1 



where new extreme points are obtained from intersecting the above planes. Next we calculate 
C-y for all the newly obtained extreme points and compare them with each other. Some easy 
calculations gives the minimum value of the parameter which is independent from A 

6 + 2(c-a) 



^'^'^^"^'"^ 9(8- 2a + 6 + c)' 
then the critical value of the parameter r is obtained as 



(9-124) 



-6 + 2(a - c) 
2 + b-c 



(9-125) 



For a — b — c — 1 the parameter r reduces to Tc = — 3 corresponding to the well known 
reduction map. On the other hand, EW (9-116) must have positive trace under any product 
state |7)(7|. Thus the introduced r in EW must satisfy 



r > rr 



-3(3 -a) > -6 + 2(a - c) 



(a + b + c-l) ~ 2 + b-c 
where the inequality is satisfied for all value of < a < 2 and a > b > c. 



(9-126) 
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10 Some separable states at the boundary of separable 
region 

Here we introduce some set of separable states as 

Pm ^^\^km){^km\ = X)|0(^l ® \l + m){l + m\, 

k I 

P'm-H \i^mk){^mk\ = E C^"*^'"'^ 10 I ® K + k){l' + k\, 

k l,l',k 

Pn-H \^nk,k){'^nk,k\ = E ^"'^'"''^10(^1 ® 1^ + k) {I' + k\, (10-127) 

k l,l',k 

where n = 0, 1,2 , m = 0, 1,2. One can show that the convex sum of po ? Po = Po' i-®' 
— jipQ + (1 — /i)pQ, is orthogonal to the optimal W\ — KWc + (1 — K)Wred-, i-e., we have 

Hence, p^^ lie at the boundary of the separable region [28]. On the other hand, one can show 
that by acting the local unitary operation Uij over VFa as {W^)ij — Uij{W^)Ulj he obtains a 
new set of optimal EW, {W\)ij, the application of which is not only to get a new set of 
bound entangled states by acting local unitary operation, but also to obtain some separable 
states {ffs)ij — UijPsUij as such which are the convex sum of separable states (10-127) at the 
boundary of separable states. 

11 Conclusion 

We have shown that finding generic Bell states diagonal entanglement witnesses (BDEW) for 
0?! (g) ci2 <8) .... (i„ systems has reduced to a linear programming problem. Since solving linear 
programming for generic case is difficult we have considered the following special cases. Also 
we have considered BDEW for multi-qubit, 2 ® iV and 3 ® 3 systems and then have considered 
optimality condition for 3 (g) 3 EW. Also, we have considered an n-d condition over 3 (8) 3 
BDEW and have obtained this condition for some special cases exactly. We have defined 
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extensive group of nd-BDEW by combining critical EW and the reduction map (each with 
special coefficients). Then we have defined the Bell decomposable bound entangled state and 
have considered detection of this state with optimal BDEW and a general BDEW. Finally, we 
have considered Choi map as an example of BDEW. Optimality and non-decomposibility of EW 
for multi-qubit and 2<S)N as well as EW for generic bipartite di <S> ^2 systems and multipartite 
(ii (g) (^2 ® • • • ® c^n ai'c under investigation. As a physical implementation of EW we know that the 
optimization of decomposition of EW to find the smallest number of measurements possible 
for local measurement on a system can be used. Therefore to make use of this implementation 
of EW for the obtained EW's is currently under investigation. 

APPENDIX A 

Minimization of the product distributions: 

In Eq.(3-5) the Bell orthonormal states for a di (8) ^2 ® ■■■ ® c^n {di < d2 < ... < dn) have 
been introduced by applying local unitary operation on IV^qq). Let us further consider a pure 
product state I7) = |q;)i|q;)2...|q;)„. Then the product distributions can be written as 

p . . (7) =|< 7 I V'. . ■ ■ (A-i) 

It easily follows that 

< P. . . (7) < ^. (A-ii) 
On the other hand, from the completeness of Bell states: 

E \Ki.,-,J(Ki2 in \-Id^^ Id, ® - «) Idr., (A-iii) 

we have Y,. . P . (7) = 1, which leads to 

E l<7|V',,,, .„>r=^^i- (A-iv) 

The above equation indicates that if we can show that for a particular choice of \a)iS, the 
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di-number of |< 7 I -0. . >P= P. . can have their maximum value equal to 4-, then 

^ III 11,12, ■■■,in I tl,t2'---'^n ^ dl ' 

the remaining ones will be zero. 

To minimize the summation C = Y^ij Qijfij for a 3 (g) 3 system, assuming that goo = 0, let 
us first suppose that \a) = so that P^^ — |. Then we find the set | {a\Uij\P) \'^— 1 for 
different possible choices of \a) and Uij: 



\a) 



\a) 



\a} — 



1 



UJ 



^0^ 



UJ 



UJ 



I \ 

UJ 



I V'oi) , 1-002) , min(Eij Qij) = Qoi + ^02 , 



1^10), 1-020), niin(Ejj Qij) = + ^20, 



I V'li) , 1-022) , min(Eij gij) = + , 



a) 



The above relations imply that C„ 
appearing in the same row. 



/ - \ 

UJ 



1 



I V'12) , |-02i) , min(Eij Qij) = + ^21 ■ 



+ 52), where qi and 52 correspond to two of qij 



APPENDIX B 

Critical entanglement witness is optimal: 

According to the References [14, 15], an EW will be optimal if for all positive operator P 
and £ > 0, the operator 

W ^{l + e)Wc-eP (B-i) 
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is not an EW. In order to prove the critical EW given in (6-75) is optimal, we first show that 

Tr{Wc\a){a\®\a*){a*\) = 0. (B-ii) 

It just suffices to check that for the product distribution P^j =< ipij\a){a\ \a*){a*\ipij >, we 
have Poo — | , Poi — P02 , Pu — P22 , P12 — -Pzi- 

Substituting Pij given above in (B-ii), it is easy to see that Tr{Wc\a){a\ ® \a*){a*\) = 0. 
Also it is straightforward to see that there exists no positive operator P with the constraint 

Tr{P\a){a\^\a*){a*\) ^0 , V\a) 

. Therefore, there exist no positive operator r to satisfy (B-i). Hence Wc, and in particular 
Wred, are optimal. 

APPENDIX C 

Simplex method for solving multi-qubit minimization problem 

We know that simplex method is an elegant way for solving linear programming problems. 
As an example we obtain the -P00...00 -^10.. .00 constraints in Eq.(4-31), thus we have two slack 
variables which are defined as 

^i = ^-2Poo...oo+2P,o...oo(l-^) , ^2 = ^-2P„ + 2P„(l-^). (B-i) 

We carry out this procedure to transform the inequality constraints (4-31) into equality 

maximize -C^ = 2(2^-1-1) (-(1 - + (1 - •^)^„«...oo + ((2'' - 1)^; - l)Pio...oo) 
subject to uji^^- 2P,„ + 2Pj„ „„(1 - ^) 

(B-n) 

^2 — 2^^^ ~ 2P^Q gg -|- 2P(,g pp(l — g^jirr) 
-Poo...oo)-Pio...oo)^l)^2 > 0. 

Now we rewrite the first equation in (B-ii) in terms of cui and 002, making use of the slack 
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variables: 



1 , ^ (l-x)a-(2"-l)x+l (1 -x) - ((2"- l)a 
2(¥^)^-^' - ^) + 2(j-l) + ^ 2(a^-l) 

(B-iii) 



where a = 1 — For < x < the coefficients cui and a;2 are both negative. Now 

from the simplex method we conclude cui — 002 — , i.e., -Poo.. .00 — -P10...00 — i- Thus the 
minimum value of C-y = |. For < x < 1, from (B-ii), we see that the coefficient P^^ ^ 

is negative, so that -Pig 00 — 0' hence -P00...00 — Therefore, we find the minimum value of 

Cj as (C-y)min 2" ' 
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Figure Captions 

Figure- 1: Feasible region for 3 (8) 3 systems for particular choice of ?oo — 0) ?io — 
Q20 = \ ~ ^ others g's being equal, i.e., when the linear programming variables are P^, 
Pio and P20 ■ 

Figure-2: Feasible region for 3 (8) 3 systems for particular choice of q^^ — 0, q^^ — x, 

q^ = \ — X and others g's being equal, i.e., when the linear programming variables are P^g, 
^10 and . 

Figure-3: Feasible region for 3 (8) 3 Choi map for particular choice of a > 6 > c, i.e., when 
the hnear programming variables are Vi Y.\=\ Puo-, ^2 Y?k=o and Vi Y.\=o Pui- 



